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Abstract 

With the advent of large scale surveys (i.e., Legacy Surveys) it is now possible to start looking beyond the 
galaxy luminosity function (LF) to more detailed statistical representations of the galaxy population, i.e., 
multivariate distributions. In this review I first summarise the current state-of-play of the B-band global and 
cluster LFs and then briefly present two promising bivariate distributions: the luminosity-surface brightness 
plane (LSP); and the colour-luminosity plane (CLP). In both planes galaxy bulges and galaxy disks form 
marginally overlapping but distinct distributions, indicating two key formation/evolutionary processes (pre- 
sumably merger and accretion). Forward progress in this subject now requires the routine application of reliable 
bulge-disk decomposition codes to allow independent investigation of these two key components. 
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1 Introduction 

For almost 30 y ears the Schechter luminosity function (LF; 
ISchechtedll979) has been the standard tool for quantify- 
ing the galaxy population 1 . The LF is loosely based on the 
Press-Sch echter formalisation for t he primordial halo dis- 
tribution llPress fc Schechterlll974l) . Moreover the LF con- 
sistently provides a good formal fit to th e observed lu- 
mino sity distribution (LD; see for example INorberg et alJ 
2002). This consistency, between the LD and LF, ap- 
pears to hold regardless of envi ronment foe Propria et al 
l2003HDriver fc De Proprisl2003l) . The only departure from 
a pure Schechter function appears to be in the central 
cores of rich clusters, where the galaxy LD is often seen 
to show a marked upturn at the giant-dwarf boundary 
(Mb w — 16 mag). Perhaps the most well known example 
is the central LD of the Coma cluster (e.g ^lTj^ntham^ggjl: 
iBeiiersbereen et aJl2002t lAndreon fc Culliandrell2002l and 
references therein) . The most plausible explanation is that 
the core contains an overdensity of giant and dwarf el- 
lipticals bolstering both the bright and faint-end of the 
core clust er LF. For example th e more extensive Coma 
survey by iMobasher et all l)2004t) recovers a flat and in- 
variant LD/LF (a = -1) to M B « -14 mag. The 
phenomen a of an upturn in the LD, has also been seen 
in Virgo Jlmpev fc BothunI Il988l : iTrentham fc Hodgkinl 



1 The Schechter function: d(<j>) = <f>* (j^) a e { ~ ~ ] d(£) has 
three key parameters, L* the characteristic luminosity where 
the exponential cutoff cuts in, <f>* , the normalisation at this 
characteristic luminosity, and a, the faint-end slope parame- 
ter. A value of a = —1 implies equal numbers of galaxies in 
magnitude intervals, a more negative (or steep) value implies 
numerous dwarf systems. 



120021) . A963 JPriver et ajlll994l) . A868 jDriver et all|2003h 
A2554 l|Smith et alll997l ) and A2218 ijPracev et aJl2004l) 

for example. However, apart from these "active" core en- 
vironments, the overall LDs from the field, to the local 
group, to the local sphere, and near & far rich clusters 
all consistently follow a smooth LF within the luminosity 
ranges probed. Fig. shows an (incomplete) summary of 
b, B, V or g-band field and cluster LFs colour corrected to 
the Johnson B filter. 

The main point to take from Fig. is that the global 
and cluster LFs each show a broad but overlapping range 
of distributions. Clearly one cannot reasonably argue for 
any significant variation between the global and overall 
cluster environment on the basis of these data. Studies 
based within the same survey data, for exa mple the two- 
degree field galaxy redshift survey study bv lCroton et alJ 
(2004), generally find fairly subtle changes with environ- 
ment. Hence it seems that the variations seen in Fig.Qindi- 
cates an unspecified systematic error in the various studies. 
The most lauded of these is the unsavoury topic of surface 
brigh tness selection effects llDisnevll97fAllmpev fc BothunI 

1997) . The concern is that the galaxy population at each 
luminosity interval occupies a range in surface brightness 
(or size). Surveys with shallow detection isophotes may 
miss both light from a galaxy's halo, as well as e ntire galax- 
ies (s e e for example | Sp ray_berrv et alJll997l and iDalcantonl 

1998) . C ross fc DrrvenT2002Tr "ex'plored this possibility in 
detail and demonstrated that indeed surface brightness se- 
lection effects can play havoc with the recovered Schechter 
function parameters and reproduce exactly the kind of 
variation seen in both the global and cluster LFs of Fig. 

More recently a number of papers have identified a 
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Figure 1. Various luminosity functions as measured for the 
global environment (upper) and cluster environment (lower). 
Data are taken from Table. 3 of Liske et al (2003), Table. 1 
of Driver & De Propris (2003), Table. 2 & 3 from Blanton 
et al (2003) and Table. 2 from Driver et al (2004). The solid 
lines show the regions over which the luminosity functions have 
been fitted and the dotted lines the extrapolations. The clus- 
ter luminosity functions have all been arbitrarily normalised to 
4>* = 0.0161 galaxies per h 3 Mpc -3 . The main point is that 
until the systematics are resolved one cannot draw any reason- 
able inference other than the field and cluster LFs are broadly 
consistent. 

clear luminosity-surface brightness (or size) 2 relation for 
field galaxies b ased on divers e datasets including: the Hub- 
ble Deep Field toriveilll999l): the t wo-degree Field Galaxy 
Redshif t Survey dCross et al.ll200ll); the Sloan Digital Sky 
Survey jBlanton et al.ll200ll : IShen et al JI2003T) : and a very 
local i nclination and dust co rrected sample of late-type 
disks dde Jong fc Lacevll2000|) . These studies consistently 
show that low surface brightness is synonymous with low 
luminosity — with a few nota ble exceptions as typified 
by Malin 1 jBothun et alll983) and the faint second disk 
surrounding NGC5084. 

To fully resolve the potential impact of surface 
brightness selection effects one must consider the joint 
luminosity-surface brightness distribution. This has been 
advocated in the past, not so much to compensate for se- 
lection bias, but to preserve the size (or surface brightness) 
informati on which may be of interest in its o wn right (see 
ICholoniew ski 1985 and lSodre fc Lahavlll993l for instance). 
This latter point is illustrated in Plate 1, where I show an 
example LF for a nearby volume limited sample and images 
of the actual galaxies contributing to the LF. Clearly much 
information is lost when one replaces these images with 

2 Luminosity, size and surface brightness are related by Mhlr = 
M + 2.5 log 10 [27ri? 2 ILR ] + 36.57 where fj, e is the effective sur- 
face brightness, M the absolute magnitude, and -Rhlr the semi- 
major axis half-light radius in kpc, hence the luminosity-surface 
brightness relation can be readily transformed to a luminosity- 
size distribution and we use the acronym LSP to indicate either. 
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Figure 2. The surface brightness distribution of galaxies at 
various luminosity intervals (as indicated). The curves show the 
Gaussian fits to the recovered joint luminosity-surface bright- 
ness distribution of Driver et al (2004). The shaded region de- 
notes the limits at which strong selection effects are likely to 
impact upon the observed distributions. Generally the distribu- 
tion is narrow and constant for the brightest galaxies (aka Free- 
man's Law) and then broadens towards lower surface brightness 
for lower luminosity systems. 

three simple numbers. It is for these reasons — the need 
to accommodate selection bias and the desire to explore ad- 
ditional parameter space — coupled with the abundance of 
data that now moves us beyond the simple LD /LF to start 
exploring multivariate distributions. Here I introduce two 
such distributions, the luminosity-surface brightness plane 
(for the reasons stated above) and the colo ur-luminosity 
plane which is also of topical interest (e.g., iBaldrv et alJ 
2003; Ho gg et al.ll2004l and references therein). 



2 Multivariate distributions 

To construct multivariate distributions requires an ex- 
tensive wide area survey (> 10°), to reasonable depth 
{Hum > 25B mag arcsec -2 ), with reasonable res- 
olution (FWHM = 1"), wavelength coverage (e.g., 
some of ub j BV grRH c z J HK 3 etc) and spectroscopic red- 
shifts/distances. The most notable catalogues for this pur- 
pose are: the two-degree Field Galaxy Re dshift Survey 
(2dF GRS; bjR F , 1800 sq deg, 250K z's; IColless et all 
1200 ID . the Sloan Digital Sky Survey (SP SS; ugriz, 10000 
sq deg, 1M z's; IStoughton et alJ 120021) . the Millennium 
Gal axy Catalogue ( MGC; B+SDSS, 37 sq deg, 10K 
z's; iLiske et all 120031) and the two Micron All Sky Sur- 
vey (2MASS: JHK , all sky, SDSS+2dF+6dF+MGC z's, 
IJarrett et all2003l) . The MGC, although the smallest in 
area, is also the deepest (fiu m — 26.05 mag arcsec -2 ), 
highest resolution (FWHM = 1.25") and most complete 

3 We specifically limit ourselves to the optical/IR regime but 
note the existence of all sky HI, X-ray and far-IR surveys 
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survey fsee lLiske et all2003llCross" et al l2004l : lDriver et all 

2004). It also overlaps with the other three surveys and 
hence provides a "best of all worlds" hybrid dataset — for 
example 50 per cent of the ~ 10, 000 MGC redshifts de- 
rive from the 2dFGRS or SDSS, extensive optical colour 
coverage from SDSS, and partial near-IR coverage from 
2MASS. The MGC 4 contains 10,061 resolved galaxies with 
12.5 < Bmgc < 20 mag with 95 per cent complete redshift 
coverage. All galaxies have been analyse d with a variety of 
softw are packag es including SExtrac tor feertin fc Arnoutsl 
Il99fj) . GIM2D dSimard et al]l2002l) and eyeball classified 
to Bmgc < 19 mag. 



2.1 The luminosity-surface brightness plane 

The luminosity-surface brightness plane (LSP) is of partic- 
ular interest because it enables one to compensate for both 
luminosity (Malmquist bias) and s urface brightness s elec- 
tion effects (aka "Disney bias"). In lDriver et all i2004i) the 
LSP is derived for the MGC, which provides the most ro- 
bust current estimate. The MGC LSP analysis used the 
joint luminosity-surf ace brightness S t ep- W ise Maximum 
Likelihood method of ISodre fc Lahavl lll993l) and incorpo- 
rate into this tracking of 5 selection boundaries relevant to 
each individual galaxy (i.e., maximum & minimum observ- 
able size & flux and minimu m observable central surface 
brightness for detection, see lDrive3ll999l) . An additional 
feature is the derivation of individual K-corrections us- 
ing the combined MGC and SDSS-DR1 colours (uBgriz). 
Fig.|2]shows the data as a series of Gaussian fits across the 
LSP at progressive intervals of absolute magnitude. The 
thicker weighted lines shows the surface brightness dis- 
tribution for the most luminous galaxies and the fainter 
lines for the dwarf regime. Two facts leap out. Firstly 
that the distributions are bounded (the Gaussian fits have 
good x 2 ' s ) broadening towards lower luminosity. Secondly 
the peak of the distribution moves towards lower surface 
brightness for lower luminosity systems. In other words low 
luminosity systems apparently show greater surface bright- 
ness diversity than giant systems. However this can also be 
interpreted in term s of the Kormendy relation for spheriods 
jKormendvlll977l) and Freeman's Law for disks dFreemanl 
These two classic studies unveiled distinct relations 
for the structural properties of spheroid and disk compo- 
nents. The Kormendy study found that the more luminous 
the spheroid the lower its central surface brightness. Con- 
versely Freeman's study found that all disks, regardless 
of luminosity, have a constant central surface brightness 
of li°B MQC = 21.65 ± 0.3 mag arcsec" 2 . The MGC result 
shown on Fig.|2]are for the combined bulge+disk systems. 
Around L* the effective surface brightness for spheriods 
and disks is fairly close — a long-time nagging coincidence. 
However moving towards lower luminosity the trends for 
spheriods and disks diverge leading to the broadening of 
the global surface brightness distribution. To investigate 
further hence requires separating out these two structural 



4 The MGC imaging and basic catalogues are available from 
http://www.eso.org/~jliske/mgc/ (additional catalogues in- 
cluding redshift, morphological and structural parameters are 
available on request from spd@mso.anu.edu.au). 
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Figure 3. The surface brightness distributions of bulge (upper) 
and disk (lower) components. The vertical dashed line shows the 
expected surface brightness for spheriods at L* (Mb ~ —19.6 
mag) and the dotted curve the Freeman distribution for disks 
systems. The shaded regions show the approximate selection 
boundaries, (lower) as above but for the disk components. 



compone nts via 2D bulge-dis k decomposition. Here we use 
GIM2D JSimard et al.l 12002) and Fig. 01 shows the data 
of Fig. 2 subdivided by structural component. The dot- 
ted line shows the original Freeman distribution which re- 
mains relevant today, albeit with a far bro ader dispersion 
than originally reported fsee lFreemarJl97oT) . It would seem 
that galaxies consist of two principle components (pre- 
sumably formed via two mechanisms: merging and accre- 
tion/collapse ?) and to unravel these two phases in detail 
must require robust bulge-disk decompositions of extensive 
samples over a variety of epochs. 
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2.2 The colour-luminosity plane 

The next most obvious key global parameter, after lumi- 
nosity and surface brightness (size) is colour, and in partic- 
ular the rest (u — r) which straddles the 4000A-break and 
hence a crude indic ator of the curren t star- formation rate. 
iBaldrv et all <2003ft and iHogg et alJ <2004l) have recently 
studied this plane extensively with SDSS data and demon- 
strate clear bimodality of the colour distribution. Fig. 21 
shows this trend for the 10k galaxies of the MGC (using 
SDSS colours) . Fig. [I] also shows this trend for the bulge 
and disks separately. To obtain the bulge colour we use 
the SDSS PSF magnitudes and to obtain the disk colour 
and we remove the bulge colour component from the global 
colour to reveal the disk colour . We now see that the bi- 
modal distribution can readily be explained in terms of 
predominantly red bulges and blue disks. This component 
segregation implies distinct stellar populations with dis- 
tinct evolutionary paths. Bulges must contain old stellar 
population and disks intermediate or young populations. 
Again this follows conventional wisdom but highlights yet 
further the important of bulge-disk decomposition and the 
need to study the component properties of galaxies rather 
than the global properties. 
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3 Future prospects 

The two distributions outlined above, both suggest that 
the well known bulge and disk components of galaxies fol- 
low distinct trends in both the surface brightness (size) 
distribution and colour distribution. This is of course not 
particularly new, however what is exciting is our ability to 
quantify these distributions and trends in detail for large 
statistical samples, and to extend this kind of structural 
analysis to higher redshift. In particular the data resolu- 
tion and signal-to-noise of the ground-based data discussed 
above is compar able to that available with the Hubble 
Space Telescope |Driver et all Il995al iDriver et al Il995bl 
IDriver et alll998bilDriveJll999n and the upcoming James 
Webb Space Telescope. There is nothing, other than hard 
diligent work, to prevent us from quantifying the evolution 
of these distributions across the entire path length of the 
universe. However three further issues are worth raising: 

1) Which wavelength is optimal for structural stud- 
ies of galaxies ? 

2) How might we push back the boundaries into 
the dwarf regime ? 

3) Can we connect structural measurements to the 
properties of the dark matter halo ? 



3.1 The near-IR 

Traditionally, almost all nearby galaxy catalogues have 
been based on flux-limited observations through an opti- 
cal B or blue bandpass filter (~ 400 — 450 nm), for exam- 
ple: the RC3; the Two-degree Field Galaxy Redshift Sur- 
vey; the Millennium Galaxy Catalogue etc. This has been 

5 i.e., {u-r) D = -2.51og[10(- ' 4uT '-10(- - 4 ( ,1 - r ) B ']-r T - 
2.5 log[10^ 1— S / T )] where the filter subscript refers to total Pet- 
rosian magnitude (T), disk (D) or bulge PSF magnitude (B). 



Figure 4. (upper) The bimodal distribution of rest-(u — r) 
colour (white) subdivided into bulge (red) or disk (blue) com- 
ponents. It is clear that the bimodal distribution is is really a 
bulge-disk dichotomy, (lower) the same data shown according to 
rest-(« — r) colour and absolute (B) magnitude. The sample is 
not volume corrected but nevertheless shows that bluer systems 
are typically of lower luminosity. 



driven by technological/commercial necessity with flux de- 
tectors typically optimised to the spectral response of the 
human eye (400 — 800nm). However the most physically 
meaningful bandpass in which to observe a galaxy is in the 
near-IR (rest-// band or 1.65pm). This is mostly because 
the stellar population that dominates the total stellar mass 
— and therefore best traces a galaxy's gravitational poten- 
tial — is the long-lived low- mass population which emits 
in the near-IR (see for example lGavazzi. Pierrini fc Boselli I 
Il99rjj . This is most clearly demonstrated by the obviously 
smoother appearance of a galaxy in the near-IR than in 
progressively bluer wavelengths (see upper panel of Fig. |S] 
showing a montage of images for M51 in a variety of fil- 
ters — the near-IR images (right most) are significantly 
smoother. The flux and shape of a galaxy in the near-IR 
is most dependent on the older relaxed stellar population 
and therefore a better tracer of the underlying potential. 
Conversely the flux and shape in the optical is linked to 
the young stars and therefore dependent on the current 
and possibly transient star-formation rate. Both optical 
and near-IR data are important if one wishes to under- 
stand galaxy formation and evolution. The near-IR how- 
ever appears to be the optimal filter for the investigation 
of the structural properties. The other great advantage is 
of course the minimisation of the impact of dust obscura- 
tion. This is illustrated in the main panel of Fig. which 



© Astronomical Society of Australia 0000 



Beyond the galaxy luminosity function 5 




0.5 1 1.5 2 2.5 

Wavelength (nm) 



Figure 5. An illustration of the advantages of the near-IR. M51 images in UBVRIJHK respectively are shown along the top. The 
main panel shows a spectrum of the night sky, a spectrum of a galaxy before and after star-burst and the location of the B and H niters. 
The extinction curve is also show. At shorter wavelengths one has to contend with the vagaries of dust and star- formation. At longer 
wavelength images are smoother, less affected by dust and star-formation. For detailed structural analysis the longer wavelengths are 
clearly optimal . 



shows the location of the B and H band filters superim- 
posed on the night sky spectrum (dotted line), a galaxy's 
continuum before and after star- burst (solid lines), and 
the dust attenuation curve (long dashed line). The im- 
pact of star-formation and dust is clearly less in the near- 
IR. The upcoming near-IR facilities and in particular the 
UKIRT / WFC AM and VISTA will have the capabilities 
to provide exactly the kind of wide, deep, high resolution 
data required for the comprehensive structural analysis of 
nearby galaxies. 

3.2 The dwarf regime 

The space density of dwarf galaxies remains elusive. 
Figs. Q & |2l shows that the MGC can only sample with 
credibility to Mb ~ —16 mag at which point both limiting 
statistics and th e high and low sur face brightness selection 
limits bite fsee iDriver et all I2004T) . Fig. HJ illustrates this 
by showing the MGC galaxies on an absolute magnitude 
versus redshift plot. The data are of course bounded by 
the B — 20 mag limit which highlights the rapidly dimin- 
ishing volume observed for low luminosity systems. One 
way to overcome this is to simply conduct ever deeper red- 
shift surveys (as indicated on the figure) . However this has 
a diminishing return as the number of galaxies one must 
observe to find one low luminosity system becomes unrea- 
sonable. One possible way forward is to use photometric 
redshifts to pre-select low-z candidates and then follow-up 
only these systems. However the accuracy of photomet- 



ric redshifts at low z is poor (although improved by near- 
IR colours, see lBolzonella. Miralles fc Pell6l2000l for exam- 
ple). To overcome the surface brightness selection limits 
(both high and low) the source data must be improved to 
probe to very high resolution (FWHM < 0.5") and very 
deep isophotes (fiB » 26 mag arcsec -2 ) over wide ar- 
eas (30+ sq deg). No such survey exists but facilities such 
as SUBARU/SUPrime and Magellan/IMACS just about 
have the capability to achieve such a survey. The alterna- 
tive method is to observe the very local galaxy population 
(i.e., the Local Sphere of Influence, defined as that within 
10 Mpc) and obtain direct distance measurements rather 
than redshifts. 

3.3 Physics of the LSP ? 

The LSP may have the potential to connect key observables 
(luminosity and size) to the fundamental underlying physi- 
cal properties of bulge and disk systems (mass and angular 
momentum). I n various studies of the formation of disk 
systems, (e.g.. jFall fc Efstathioul Il980t iDalcanton et. alJ 
Il997l iMo et alJ Il998f) the dimensionless spin parameter 

(A = J\E^\G- X MZJ.. iPeeblesI ll96Sh is directly related 
to the scale-length of the disk. The spin parameter re- 
flects how close the halo is to a rotationally supported 
system and is a key p arameter monitored by t he nu - 
merical simulations (see ISteinmetz fc B artelmarml (0,995); 
ICole fc Lacevlll996HVitvitska et alJl2002t Mailer. Dekel & 
Sommerville 2002 for example). The pivotal idea (here 
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Figure 6. The difficulty of quantifying the faint-end of the lu- 
minosity function is highlighted in an M v z plot such as this 
one. The data points are from the MGC which is limited at 
B = 20 mag. The yellow box marks its boundary and it barely 
contains any volume for very low luminosity systems. One pos- 
sible way forward is to simply push progressively deeper with 
the upcoming multiplex spectrographs (e.g., AAf! on the Anglo 
Australian Telescope or GMOS/KAOS on Gemini). The main 
problem with this approach is the sheer numbers of objects. To 
circumvent this one can envisage implementing a photometric 
redshift cut first to pre-select candidate low-z objects. The left 
side indicates the distribution of galaxies from the local group, 
Fornax and the Milky Way globular cluster distribution as in- 
dicated. . 

echoing the toy model of Ide Jong fc Lacevl |20Q(J) starts 
with the premise that the baryons are coupled to the dark 
matter halo, because of this the luminosity (generated by 
the baryons in the form of stars) can be related to the sys- 
temic mass and the rotation of the stars/gas can be related 
to the systemic angular momentum. Given this premise, 
which is intimated by the Tully-Fisher relation, one can 
analytically relate A to luminosity and surface bright- 
ness (or size): A oc T, ef 2 f L~i + ^ (from Ide Jong fc Lacevl 
2000), where £ e // is the effective surface brightness, L 
is the intrinsic luminosity in some filter and 7 is the de- 
pendence of luminosity on th e mass-to-light ratio (equal 
to 0. 69 in B or 1.00 in H iGavazzi. Pierrini fc Boselli I 
1996). Numerical simulations consistently find that the 
distribution of the spin parameter is a log Normal dis- 
tribution which is globally pr eserved through hier archi- 
cal merging (see for example IVitvitska et alJ I2002T) this 
yields: £ e // = L°g 5i or fJ, e ff = 0.54Ms. Hence the gra- 
dient of any luminosity-surface brightness relation bears 
upon the relation between mass and light and the dis- 



Figure 7. A summary of available LSP data drawn from a 
variety of sources. The red line marks the credibly mapped area 
and the cyan line shows the expectation from de Jong &; Lacey 
(2000). This appears to follow the data remarkably well. 

persion upon the breadth of the spin distribution. Fig. [7| 
show s the _B-band LSP fo r a variety of samples as indicated 
(LG lMateclll99Sl : HDF. lDriveilll999l : MGG. iDriver et all 
120041: MW G Cs, Harris et al (priv. comm ) ; Loc al Sphere 
of Influence. [JjsrjeiTjJ^inj^ l2000l : LSBGs, 

Ide Blok. van der Hulst fc Bothunl Il995l) . The solid lines 
show the approximate expectation as argued above and 
show remarkable agreement with the data - in detail the 
observed size distrib ution is marg i nally n arrower than sim- 
ulations predict (see lDriver et all ll2004) ). It is also worth 
noting that systems which form via merging (i.e., bulges) 
and via accretion (i.e., disks) are also predicted to show dis- 
tinct A distribution s (see for example I Vitvitska et al.l2002l : 
iMaller et aiT 2002). At the moment far more data and de- 
tailed simulations are required however this connection is 
clearly promising and could ultimately result in a galaxy 
equivalent to the Hertzsprung-Russell diagram, allowing a 
meeting ground between numerical simulations and survey 
observations. 



4 Summary 

The galaxy luminosity function is not the only fruit, and 
with the many legacy datasets becoming available the time 
is now ripe to move beyond the LF and explore multivari- 
ate distributions. Here I've presented two: the Luminosity 
surface brightness (size) plane (LSP); and the colour lu- 
minosity plane (CLP). Both planes show that disks and 
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bulges form distinct but overlapping distributions presum- 
ably indicating secular evolution of these components, i.e., 
two mechanisms and two timescales. This finding argues 
for the community to move away from global measure- 
ments and start to measure the properties of these distinct 
components independently. We argue that this is best done 
in the near-IR and should be a key focus of upcoming IR 
facilities such as VISTA (low-z) & JWST (high-z). Per- 
haps most important of all the LSP appears to provide a 
direct meeting ground to the numerical simulations. This 
last point is by far the most important as its from the cross- 
talk between simulations and observations that real insight 
into the processes of galaxy evolution and formation will 
come. 
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Figure 8. PLATE1: The global Galaxy Luminosity Function (red line) condenses the available information of galaxies (images) into 
three crucial numbers: the characteristc luminosity (Af*); the absolute normalisation (</>*); and the faint-end slope (a). Although the 
schechter parameterisation is more often that not a remarkably good fit, one cannot help but feel that too much important infomation 
may have been lost, for instance the sizes and bulge-to-total parameters. 
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